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Abstract . A brief critique of spherically symmetric conventional 
modulation theory is supplied. Estimates are made of the cosmic 
ray intensity at high solar latitudes. Direct evidence for 
significant off-ecliptic cosmic ray gradients is reviekred in 
support of the requirement for an off-ecliptic spacecraft mission. 
The possibility of measi'ring the galactic spectrum is discussed. 


Spherically symmetric modulation theory and its problems . 

Cosmic ray moJuiation arises because motion of the energetic particles 
along the spiral interplanetary magnetic field lines is controlled by 
scattering due to magnetic i rreguJari ties which are being convected 
outwards by the solar wind. Inward diffusion is balanced both by outward 
convection and energy loss of the particles as they suffer adiabatic 
deceleration in the expanding wind. A Fokker-Planck equation expressing 
these effects in a spherically symmetric, steady-state situation is 
(Parker 1965, Gleeson and Axford I 967 ) 

7 ^ If 1 - ! ; ^ '“™> ") 

Here U is the differential particle number density at position r and 

T + 2T 

kinetic energy T,a= j --' with as rest mass energy, V is the solar 

I + I o 

wind velocity and K,-, the effective radial diffusion coefficient, is given 
by Kf ■= Kji cos^^ + Kj^sin^^ for K jj ar.d K , respectively the parallel and 
perpendicular diffusion coefficients with •= cos'* f l®.| l£.l )• 

Particle streaming ^ is due to diffusion in the wind frame plus an 
additional term involving a Lorentz transformation to the rest frame via 
the Con^ton-Getting factor C such that 
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S “ cyu - K.grad U 


( 2 ) 


with ^ ^ ^ 3 y ^ (oTll) and K as the tensor diffusion coefficient. 

At inagnetic rigidities exceeding about 1 GV, radial streaming is 
usually negligible and under these circumstances the Fokker-Plancfc 
(1) can be integrated to give a relati<»t between the Galactic density 
U(rg, Tg) and the observed density U(r,T), 

U(r, T) = U(rg, Tg) exp (- ^ dr) (3) 

for spherical symmetry with Tg determined by the effective outer bound 
to the interplanetary scattering process. 

It is important for Astrophysics to know the energy spectra of 
various nuclear species of cosmic rays in the galaxy and a common method 
for achieving this "demodulation" is as follows: 

(a) Estimate the Galactic electron spectrum from radio synchrotrtxi data. 

(b) Compare the near-earth electron spectrun with the Galactic spectrum 
to find the magnitude and rigidity dependence of /r*^® (^/Kj.).dr. 

(c) Use the results of (b) to correct the near-earth proton and heavy 
nuclei spectra for modulation. 

Various objections can be raised to this scheme. First, the average 
Galactic electron spectrum derived from radio data may not represent the 
local spectrum and in any case one must be sure of the model enabling the 
effects of local, cold, absorbing interstellar clouds (Goldstein et al. 
1970a) to be taken into account. Second, there is no adequate theoretical 
explanation as yet for the magnitude, rigidity dependence and radial 
dependence of K^, as witnessed by investigations of solar proton diffusion 
(e.g. Webb et al. 1973) end the measurement of cosmic ray radial gradients 
(<l0i/AU) which are much less than those expected on the basis of 
theoretical values (e.g. Webber and Lezniak, 1973) • It is important to 
know the r dependence of the integrand in (3) since adiabatic deceleration 
is inversely proportional to r. This importance becomes apparent in our 
third point, based on the work of Goldstein et al. (1970b) and Gleeson and 



Urch ( 1971 ) for certain choices of r dependence and reasonable 

magnitudes of find that adiabatic deceleration may completely 
exclude all <100 NeV/nucleon particles from the inner solar system. 

Thus no certain knowledge of the lowest energy Galactic primaries is 
possible. Fourth, the cause of time variations in the modulation is 
unknown, with not enough solar cycle variation in the solar wind 
par^Kters being observed near earth to account for the known changes in 
the cosmic ray flux (e.g. Hedgecock et al. 1972). For example, the 
power spectrum of magnetic irregularities at 10'*^ Hz should change by a 
factor >2 in order to account for the observed modulation of 1 GV particles 
between I 965 and 1969 according' to rescmant, wave/particle interaction 
theory while in practice the change is <I0« (Hedgecock 1975)- Off- 
ecliptic control of modulation via the effects far beyond 1 AU of solar 
streams emerging from the zones of maxiimjm solar activity may be the 
only way to explain the cosmic ray 11 -year cycle. Hence careful study 
by off-ecliptic spacecraft is required. 


The off-ecliptic route to the boundary. 


It has been thought that a direct determination of the Galactic cosmic 
ray charge and energy spectrum can be achieved employing ecliptic plane 
spacecraft trajectories to the outer planets. In this way the problems 
of section 1 are all by-passed. However, the boundary to modulation 
could be as far as 100 AU and the low, measured cosmic ray density 
gradients seen by Pioneer 10 render this approach uncertain. Alternatively, 
we note that cosmic rays have an easier inward motion over the solar poles 
where the geometric path length along the interplanetary field lines is 
much shorter than in the tightly wound spiral regime of the equatorial 
plane (Lietti and Quenby I 968 ). An appropriate Fokker-PI anck equation 
for the steady state which relinquishes the requirement of spherical 
symmetry and takes into account the spiral geometry is then 

if K j ^ 0 for path length s along a magnetic flux tube of area A, 

Allowing = K^r sin^Q where 0 is solar latitude yields 
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and if q o and ^/Ug Ve at 1 GV, we obtain the following table for 
the percentage residual modulation at 1AU: 


0 = 90® 0 = 30® 

“/Ug (at 1 GV) 63i 22 % 


0 -= 0 ® 


9 % 


Thus on a single model for scattering, a spacecraft passing over the poles 
at 1 AU may see '^90% of the unmodulated intensity and therefore get a better 
measurement of galactic conditions than is available at Jovian distances. 


3. Direct evidence for off-ecliptic gradients 

Observations confined to the ecliptic plane can only reveal the existence 
of off-eciiptic effects by noting modifications to cosmic ray streaming 
from the expectations of the spherically symmetric model case. Equation 
(2) in a more general form is 

i •'i, <17 -i:< it' <« 

where v is energetic particle velocity and u is cyclotron frequency. It 
has been assumed that direct slippage of particles across field lines makes 
only a small streaming contribution ^ ^ being time to 

travel one parallel mean free path). Furthermore, we assume the short- 
circuiting by scattering of that perpendicular gradient which would cancel 
out the anisotropy due to (£x^)^drift in the non-scattering limit when 
Liouvi lie's theorem applies to the particle intensity. At high rigidities, 
some few to a hundred GV, the radial streaming is negligible over long 
periods with the third term on the right of (b) cancelling out on average. 
Then the first two terms combine to give the streaming from the east or 
1800 hr LT anisotropy. When, however, the anisotropy is studied in 
practice as a function of sign of the interplanetary field sector structure, 
two effects of the third term become apparent. A north-south anisotropy 

X II 

arises due to ( -r— ) , x B, or the effect of the radial gradient and an 

ecliptic plane anisotropy arises due to ^ effect of off- 
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ecliptic gradient. Hashim and Sercovitch (1972) find G 2 “ S.5 R~^*^%/AU 
directed north -»• south in 1^7/68 for the tatter effect, possibly 
physically resulting from the eicoess northern hemisphere solar activity 
suggested at that time by coronal green line S303A emission. 


The previous discussion refers to the first derivative of density, 
but studies of the second harmonic of the cosmic ray intensity can reveal 
the presence of a rising or falling, symmetric, off-ecliptic gradient via a 
dependence on the second derivative. Lietti and Uuenby (1988) essentially 
use a version of (5) for the rising gradient case to predict a second 
harmonic^wi th direction of maximum perpendicular to ^ and amplitude 
a2 “ J S at rigidity P, cyclotron radius p. This 

expression is in reasonable accord with ooservations although Nagashima 
et at. ( 1972 ) claim that a cylindrical pitch angle particle distribution 
about S[ with a different physical cause better fits cosmic ray anisotropy 
data. 


With finite K , the symmetrical gradient will either feed particles 
into the equatorial plane or draw them off to higher latitudes, thus 
setting up radial streaming. A correction to the Fokker-Planck (1) is 
employed by adding „ 

dlv ( 7 ^ ( ) 1 to the right hand side. Dyer 

u C _ 

et al. (1974) in particular evoke this streaming for a falling gradient to 
explain a sunward flow at > I GV seen by a satellite detector which is 
too large to be explained by any energy loss effects in a spherically 
symmetric model. These authors require maximum modulation over the 
sunspot zones with maridianal flow patterns set up to draw particles down 
from 0 *“0® and up from 6 =90®. Cecchini et al. (1975) have developed 
a computational model to confirm the above model. Chief features are: 


K = KqBP f(0) I 


"'o 

(Vr„)’ 

(Vr„) 


r < 1 AU 
r > 1 AU 

r < 1 AU 
r > 1 AU 


f(8) = 1 + cos 8 (-3 + 5 cos^ 0) ; 




2.2 10^^ cm2 


/sec Mev 


at 1AU ; 



e “ 0.5t tAU; V„ — const.; Pq 


■=■ 100 MV; Ug o(T + T0)"2-75 

Thus ‘=’0.05 at 1 GV. The results of employing an alternating gradient 

technique to solve the corrected Fokker-Planck, with finite 1^, is to predict 
an inward streaming 'x>0.3% in amplitude between 2 and 10~^ 6eV, a radial 
gradient ''>10%/AU at 1.1 GeV between 1 and 10 AU and a ratio 
at 1 AU varying from 2.5 at 9 •“ 0 to 0.7 at 8 = 60° for 1.1 GeV protons? 

Hence it is possible to explain the radial streaming with reasonable gradients 
and Ki/,. ratios. 


Conclusions. 


We have shown that study of cosmic ray modulation by integrating the transport 
equation outwards in the ecliptic plane, assuming spherical symmetry, 
encounters various problems. The transport processes and boundary conditions 
are insufficiently well understood, modulation may be controlled by off- 
ecliptic gradients and assymetries can have noticeable effects on solar 
equatorial plane observations. Three-dimensional study of the solar 
cavity cosmic ray distrijution is required to: 

(a) Measure off-ecliptic gradients and streaming. 

(b) Enhance understanding of the solar control of intensity time 
variations. 

(c) Gain better knowledge of boundary conditions, especially the 
possibility of measuring a near-Galact ic energy spectrum over the 
solar poles. 

Objectives (a) and (c) are satisfied by a Jovian swing-by mission 
but (d, requires a direct injection at 1 AU spacecraft for detailed time 
variation studies on solar wind and solar parameters. 


RFPnODUCIBiLlTY OP THE 
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Abstract 


Based on a modified WKB analysis of the interplanetaary irregularity spectra, 
a discussion of the radial dependence of the radial cosmic-ray diffusion coef- 
ficient at polar heliographic latitudes is presented. At 1 AU radial distance 
the parameters are taken to equal those observed in the ecliptic. In the sense 
of. a present best estimate it is argued that relativistic nuclei should have signi- 
ficantly easier access to 1 AU at the pole than in the ecii.iJtic. The reverse 
may very well be true for the direct access of very law rigidity particles. 




uotb-ttv of^ 
VAGE is poor 



I, Introduction 


The access of galactic cosmic rays to the inner solar system is regulated 
by the cumulative effec*" of the irregular fields which scatter the par- 
ticles on their way in. To determine the degree of this access . . helio- 
graphic latitudes that are significantly different from the solar equa- 
torial plane, one has to rely- on theoretical estimates. These conventional- 
ly assisne a diffusive propagation scheme. Of primary importance is then 
the spatial diffusion coefficient or, more generally, the diffusion ten- 
sor, and in particular its spatial variation in the solar system. 

In the earlier work of V61k et al . (1974) the irregular fields were assumed to 
be Alfv6n waves of solar origin, convected outwards by the solar wind. 

In a (on the average) stationary interplanetary medium with axial symme- 
try around the solar rotation axis, this led to a radial direction for 
almost all wave normals beyond about 1 AU. Assuming the effective average 
magnetic field to be in the ideal soiral field direction and the wave am- 
plitudes to vary according to the (WKB) approximation of geometrical 
optics, a radial dependence of the coeff>.c.ent for diffusion along the 
average field was calculated. As a sensitive function of the angle be- 
tween wave normals and average field, the value of this diffusion coef- 
ficient, or equivalently , of the mean free path, varies from a minimum 
at O degrees to Infinity at 90 degrees. Since in the solar equatorial 
plane the angle between the radial and the spiral field direction is 
large already at 1 A'.', and increases with radial distance, deviations of 
the wave normals from the radial direction have been discussed by Richte r 
(1974) in the context of the solar wind stream structure; see also Hollweg 
(1975) . 
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Ho do not intend to evaluate here the effects of deviations from the 
siiqple picture of V61k et al . (1974) at low latitudes and refer to a 
future publication ( Morfill et al«» in preparation) . Ka shall rather 
concentrate on the situation at polar heliographic latitudes. There the 
assumption of all wave vectors Ic being radial and parallel to the average 
field ^ leads to the minimum value of the radial diffusion coeffi- 
cient K as far as its dependence on the angle between ^ and B is 
concerned. Thus* within the HKB approximation for the radial development 
of the scattering centers, this provides a lower limit for . Such a 
lower limit is of interest, if- it can be argued, or at least be specu- 
la t d in a reasonable way that even in this case there is essentially un- 
impeded access to eibout 1 AU for a significantly larger part of the galac- 
tic energy spectrum than in the solar equatorial plane. 

He shall give a short discussion here, motivated by three considerations. 

The first one is that in two earlier publications this author was associated 
with ( V61k et al . , 1974; V61k , 1975), a rather different result was be- 
lieved to hold true. Another reason is given by the present discussion of 
an ex-ecliptical probe to explore the sun and the interplanetary medium. 

The final reason is that via easy access at the poles much larger regions 
of interplanetary space might possibly be populated by particles of galac- 
tic origin. 

In the next section we shall present che general behavior of k. in a modified 
HKB analysis, using simple approximations to two rather different measured 
interplanetary power spectra. In the last section the relation to the expected 
actual situation at the heliographic pole is discussed. 
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2. WKB Analysis 


Consider the singles t, axlsyometrlCf interplanetary nedium# where all quanti- 
ties depend on heliocentric distance r and where only the ideal spiral field 
B with components 


( 1 ) 




i 


VV 

r 


yiH^‘Sin9 

V ’ 


and the power spectrum P may in addition depend on heliographic la- 
titude 9 . Here r is a radial reference level/ B is independent of 9 ; 

^ denotes heliographic long itude,^gfi^ 2.65 x 10~^ cps is the angular fre- 
quency of the son's rotation/ V is the solar wind speed (assumed to be radial 
and constant) / and f is ttie wave frequency seen by an observer at test re- 
lative to the sun's center. Then it is simple to show ( V61k et al . » 1974) 
that for Q — r ^ the radial diffusion coefficient is given by 




1 



In equation (2) we have ^ ratio of particle velocity w to the 

velocity of light c; R denotes particle rigidity; v^ “ B • (47C^ )"^ is 
the (vectorial) Alfvfen - speed, where ^ is the average solar wind mass den- 
sity;^ - is the cosine of the particle's pitch angle, where w^^ is 

the velocity component parallel to B; r^ ■ 1 AU. The ^^1 i f ic^t^t^n ^^ ^c^^ 

ouiolnal page is ^OOR 
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for th« Al£v6n wave amplitudes Sb^ 


ISSjff (r) I 

( 3 ) 2 - 


/V't Vav (^) \ ^AiT M 

\ V+ (r) / 'V W 



is independent of wavevector and b' • If both and r are larcje compared 

to the Alfv&nlc critical radius (S$ 20 in the equatorial plane) , then 

the approximation on the far r.h.s. of equation (3} holds in this lowest 

order WKB approximation . Kov;uver, following equations (1) and (3) we have 

(at® » 90°); (r ,6)^ / B^{r,©) ow r, where is the total 

2 2 

power in the fluctuations. Thus, the possibility arises that B ‘ > 1 

beyond some radial distance in which case we expect wave propagation to be 
rather drastically altered by nonllneai. effects that ultimately should lead 
to dissipation of wave energy. To take this possible effect into account, we 
modify equation (3) so that for all r 

(3a) ^ Sb (r ) ^ B (*<) 


A similar device has been used by Hollweg (1973) and appears as a simple if 
crude way to take the inadequacy of the linearised WKB approx iaiation into 
account. 




It is clear physically that equations (1), (2), (3) and (3a) alsc 

at given r ,^ro' 5 different from their values at Os^c, ar, long as 

their dependence on Q is weak.enouah, such that for example — !• 

V dG 

For the rest of this section, however, we will consider (rj), V, ^ (rj) 
and P(f,©,rj) as given by their values in the equatorial piano. 



Observed power ^ectra at Qt&O and rs^r^ (e.g. Joktpii and Coleman , 1968; 
Bercovttch, 1971), generally approximate rather well a power law for P(f,r|) 
at high frequencies, tdille flattening at low frequencies k Qualitatively, an 
analytical form 

^nvf.) 

witti parameters C^, f^^and q <0|ls not an unreasonable representation. With 
equation (4) the general character of 1C (r) Ceui be inferred. For small enough 
rigidity R, so that . v/R> ^ we have K whereas for 

jal* v/R^ 27Cf^. The exponent of r is additively increased by unity when- 
ever the supplementary restriction (3a) comes into force. For our present 

3 2 

discussion we choose q = - /2. Thus, /C <55 co.ist (or *^r) for small R . r , 

3 4 2 

whereas K *vr (or *vr ) for large R • r . At large r, AC will 

4 4 

be »^r (or /s»r ) for any R. At fixed values of q and R, the transition of 
3 4 

Kt to the r (or r ) dependence increases with decreasing f^cf . equation (2) . 

As an aside we should mention that for this value of q and for the present 
case of k H B the quasilinear expression for K used in writing in equation 
. (2) is numerically quite satisfactory and the modifications in the region 
around |l = 0 (e.g, Jones et al . , 1973) therefore not essential. 

Numerical results are shown in Figures 1 and 2. Figure 1 is calculated using 
the spectrum of Jokipii and Coleman (1960) where, in reference to equation 
(4), Cj RJ 16 X 10 ^ • ^^(Hz) ^ ^^<5^7 x 10 Hz^’ and Alfven wave 

propagation was started at r^ = 20 Rq (solar radii) . Figure 2 uses the spec- 

_3 

trum published by Bercovitch (1971) which we roughly approximate by C «6 x 10 . 

2 1 -6 3 * 

•{Hz)"*"‘J; f ^ {5^ 7 X 10 "Hz; qsJ- /2, Although the results are not very much 
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different for both cases near cel AOt the very different f^-values lead to 

3 4 

strong differences in the «iset distance of the r (or rather r ) law. Ms 
diould mention here that %>e used here tte coovonent of the magnetic spectral 
tensor perpendicular to the ecliptic to represent the total power per fre- 
quency interval. For true axisymmetry of the spectrum all 1C values in Figures 
1 and2 should be multiplied by a factor 1/2. This is an unavoidable uncertain- 
ty 

The interesting aspect of these results is that they ioq>ly little modulation for 
relativistic nuclei, where adiabatic deceleration is small^ considering the 
value of f^ in Figure 2 as a rather extr^ae lower limit to the actual situation. 
For 1 GeV protons, for example, a lO-^ percent modulation is estimated in 
the diffusion convection approximation. 

3. Discussion 

The above results were obtained by fitting^ at l.AU^the average solar wind 
parameters as well as power spectra by the corresponding quantities observed 
at 1 AU in the ecliptic; the spatial dependence of the spectra assumed a modi- 
fied WKB approximation for the irregularities. In reality, the polar region of 
the corona may well be a large, stationary coronal hole, resulting in a some- 
what (perhaps fifty percent) larger flow speed and, possibly, somewhat more 
power in the frequency region f f^. All this would lead to a moderately in- 
creased modulation. 

On the other hand, it may very well be that the power 

in frequencies ^ is much smaller at TC/2 than near 0 = O. We have taken 
these fluctuations also to be Alfvenic which leads to the amplitude variation 
given in equation (3) .In reality, the part in the spectrum with f^ may well 
be due to the solar wind stream structure ( Goldstein and Siscoe , 1972) . If 
the latter is assumed to be absent at © = 90^, also the power at f^ would 
be absent with a corresponding decrease in modulation. In this light, also 



tbe possibility of increased scattering at larger distances due to local pro- 


duction of waves - a situation that is quite likely at 0^0 - aj^ears 
rather weak. Irregularities produced by enhanced (coi^ared to the ecliptic) 
tnrtmlence due to rcuiially increasing departures frcHs thermal equilibrium at 
tte poles should be assumed to have small scales, irrelevant for coaalc rays 
even in the case of the Firehose instability. It should be added that in con- 
trast to a pc^nilar feeling this result for fC and the consequent argument for 
modulation has little to do with the shorter gecmetrical path along ^ of a 
galactic particle to, say, 1 hU, but rather to the strong decrease with r of the 
magnetic field at the poles. 

Ihus, barring unknown new effects, the present best estimate is that cosmic 
ray access at the poles should be significantly better than at rela- 

tivistic nuclei. For very low-rigidity particles on the other hand, the sharp 
increase of tC with r occurs only at such a large radial distance that their 
direct access may be at least as strongly prohibited as in the equatorial plane. 
BOwever, for this last kind of statement^, the present estimate is not well 
suited. 
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Figure Captions 

Figure 1: The radial diffusion coefficient /C at heliographic latitude 

Q B 90 ° as a function of radial distance r in solar radii 

for various proton energies. Ihe power spectrum P(f)=C^ • f*^ 

with the values = 16 x lO ^ (Hz) ^ f^= 7 x 10 ^ Hz 

(Jokipii and Coleman , 1968) . Wave normals are assumed to be 

radial. The calculation was started at r = 20 R ^ , 

o V 

Figure 2: The same as Figure 1 with values of and f^ adapted to the spec- 

tra of Bercovitch (1971) . 
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